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Abstract 

Density functional theory and ab initio molecular dynamics simulations are applied to investigate 
the initial steps of ethylene carbonate (EC) decomposition on spinel Lio.6Mn204 (100) surfaces. EC 
is a key component of the electrolyte used in lithium ion batteries. We predict an slightly exothermic 
EC bond breaking event on this oxide facet, which facilitates subsequent EC oxidation and proton 
transfer to the oxide surface. Both the proton and the partially decomposed EC fragment weaken 
the Mn-0 ionic bonding network. Implications for interfacial film made of decomposed electrolyte 
on cathode surfaces, and LL r Mn204 dissolution during power cycling, are discussed, keywords: 
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I. INTRODUCTION 



Lithium ion batteries (LIB) featuring transition metal oxide cathodes and organic solvent- 
based electrolytes are currently the energy storage devices used to power electric vehicles. 
The high operational voltages of LIBs allow them to store significant amount of energy How- 
ever, these voltages can also induce electrochemical side reactions at electrode-electrolyte 
interfaces.- - - These lead to irreversible capacity loss, power fade, durability, and safety issues. 
Understanding and controlling these interfacial reactions are of great interest for improving 
electric vehicles. 

In this work, we focus on the basic science of one such reaction, and apply first principles 
computational techniques to study the first steps of ethylene carbonate (EC) decomposition 
on the (100) surface of spinel Li r Mn 2 04 used as the cathode material in many lithium ion 
batteries.- - - EC (Fig. [1]) is an indispensible electrolyte component for batteries that rely 
on graphitic carbon anodes. At the negative voltages needed to charge up a battery and 
intercalate Li + into graphite, EC-containing electrolyte decomposes to form stable solid- 
electrolyte interphase (SEI) films that prevent further electron leakage to the electrolyte. -~- 
As such, substantial theoretical work has been devoted to excess electron-induced EC chem- 
ical reactions in the vicinity of the anode.- - — One may argue that its ubiquity and critical 
role make EC in LIB the battery equivalent of the H 2 molecule in biology, geochemistry, 
and many solid-liquid interfacial science disciplines. The present work is indeed modeled 
after theoretical studies on water-on-mineral surfaces.— - — 

Oxidative reactions of organic electrolytes on cathode surfaces and formation of surface 
film from electrolyte decomposition products there^ - — are arguably less well understood 
than reduction on anode surfaces, especially at the atomic lengthscale. Electrolyte de- 
composition products such as LiF,— acetone,— aldehydes, 31 carbon dioxide,— i2L2i organic 
radicals,— and unidentified polymer, polyether, and carboxylic acid species^ 1 ^ - — >2a have 
been detected either in gas form or on the surfaces of Li x Mn 2 04, as well as on Li x Co0 2 and 
noble metal electrodes set at voltages comparable to those of Li x Mn 2 04. The salt used in the 
electrolyte strongly affects the surface film composition and properties.— ^IMn^^L^l For ex- 
ample, LiPF 6 yields fluorine compounds that slow down Li + transport into the cathode.— >22 

The electrode is also degraded by side reactions. After cycling power, Lf r Mn 2 04, ex- 
hibits surface morphology changes and loss of Mn and oxygen [ ons ±A^jMjM.^l>M-M.^i^M^. 
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It has been widely accepted that Mn(III) ions on spinel electrode surfaces disproportionate 
into Mn(IV) and Mn(II),— the latter of which dissolves into the electrolyte. The presence 
of acid has been shown to facilitate dissolution and cathode degradation.- In organic sol- 
vents, acid has been speculated to come from reactions between PFg~ and trace water,— 
or from organic solvent decomposition products themselves. 32 Counterions other than PFg 
cause slower Mn204 dissolution.—^ There is also evidence that the emergence of proton- 
intercalated surface regions^ or other surface phases,— not just the loss of Mn active mass 
from the electrode, is responsible for cathode capacity loss. Dissolved Mn ions have been 
shown to diffuse to the anode, degrading the SEI there and weakening anode passivation.— - — 
Thus anode and cathode degradation are closely related. Despite these key insights from 
experimental efforts, the free energy changes, activation barriers, specific surface configura- 
tions, and reaction pathways associated with Mn(II) dissolution have not been elucidated 
in organic electrolytes. 

On the computational side, studies have focused either on cathode oxide surfaces without 
solvent molecules,— - — or on anions^ 1 ^ or EC + /PC + (with an electron removed) in the 
absence of the oxide.— - — In studies of oxide- vacuum interfaces,— - — it has been shown that 
under-coordinated Mn ions on Mn 2 C>4 surfaces adopt low oxidation states, and these Mn 
have been argued to be susceptible to dissolution.— However, dissolution likely require that 
solvent molecules first directly coordinate to surface Mn ions,— which restores part of the 
octahedral coordination environment. When this happens, it is not obvious that the surface 
Mn(II) and Mn(III) ions do not revert back to Mn(IV). Explicit inclusion of solvent molecules 
on electrode surfaces would be extremely useful to investigate conjectures about surface Mn 
charge states.— 

The oxidation potentials of EC (and the similar PC)^ and electrolyte anions^^i such 
as PFg have been computed using DFT methods. These calculations assume a bulk liq- 
uid environment described by a dielectric constant e , and an intact EC molecule (no ring 
opening reactions during oxidation). The predicted oxidation potentials are at least 1.5 V 
above the maximum voltage of 4.3 V applied in many electrolyte oxidation experiments 
for undoped spinel Lia.Mn 2 04,— i^Lju^iil - n j-^^ agreement with tabulated data obtained 
using inert electrodes.- Recently, intriguing theoretical studies have predicted that PFg or 
ClOj in the vicinity of PC, with e reduced to reflect "near-electrode" conditions, lowers 
the PC oxidation potential and leads to novel low-barrier reaction pathways.— 1 ^ The anions 



3 



act as proton shuttles and/or provide fluorine atoms to react with PC. The oxidation steps 
are accompanied by bond-breaking events, which occur spontaneously during geometry op- 
timization of (say) the [EC/PFg ] + pair. However, the predicted redox potential of this pair 
(4.94 V) remains higher than 4.3 V even after the local e a is set to unity— This suggests 
that the electrode surface plays a role in catalyzing electrolyte decomposition. Despite the 
assumption of the proximity of the cathode, the electrode itself is not explicitly included in 
the calculations. Nevertheless, the predicted reaction mechanisms successfully reproduces 
organic products observed in experiments-^ 

In the present work, the effect of Lio.6Mn204 electrodes on EC oxidative decomposition 
is explicitly considered at the atomic lengthscale. The x = 0.6 Li content is roughly at the 
halfway point of charging or discharging. Our focus is on the initial reactions that modify 
and coat the oxide surface, not on the final distribution of gas/organic products. Our calcu- 
lations do not extend to the high voltage region which seems to be required for CO2, one of 
the proposed final products.— Given the inevitable mixture of crystal facet/termination in 
battery cathode materials and the myriad effects of applied voltage, state-of-charge (lithium 
content), temperature, electrode cracking, sweep rate, presence of conductive carbon black, 
dissolution of surface films formed by electrolyte decomposition, and interference from salts, 
contaminants, and even processes at the anode, adopting all experimental conditions to 
these theoretical calculations is extremely difficult. Instead, we take a basic science starting 
point and focus on EC breakdown on defect-free Li 06 Mn 2 O 4 (100) surfaces. Our predic- 
tions can potentially be verified with experiments conducted under ideal conditions, but the 
mechanism predicted and insights obtained may still be relevant to battery conditions and 
may motivate new experiments and interpretations of battery degradation processes. 

Our approach resembles the intensive study of H 2 adsorption on rutile (Oil),— anatase 
(101),— »^ anatase (001),— and many other surface facets of technologically important tita- 
nium oxide over more than a decade. As in the case of water-oxide interfaces, both ultra-high 
vacuum (UHV) and liquid-immersed electrode conditions are considered in our interfacial 
models. These conditions yield qualitatively similar predictions. Like H 2 on some Ti0 2 
facets, EC is found to chemisorb on Li r Mn 2 04 (100), breaking internal covalent bonds in the 
process. This initial decomposition step enables the subsequent oxidation reaction which 
involves electron and proton transfer to the oxide surface. This suggests that organic solvent 
decomposition and Mn 2 04 dissolution may be related. No assumption of a local e a is needed 
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in our calculations. 

The (100) surface has been the focus of most theoretical studies on spinel L^M^O^— i5£i£I 
Although the (111) surface is most prominent in LiMn204 particles,— (100) facets remain 
even after annealing at 800 °C for 100 hours.— The (100) surface has also been synthesized 
directly on strontium titanate substrates.— 

This work is organized as follows. Section 2 describes the theoretical method used. Sec- 
tion 3 discusses EC decomposition on Li 06 Mn 2 O 4 . Further discussions in Sec. 4 put the 
results on ideal model surfaces in the context of battery operating conditions and provide 
more comparison with the theoretical literature. Sec. 5 recapitulates the key findings, and 
two appendices describe the charge states of the Mn ions and the reactions of H 2 contam- 
inants and hydroxyl groups on Lia.Mn 2 04 (100) surfaces. 

II. METHODS 

Static (T=0 K) DFT+U and finite temperature ab initio molecular dynamics (AIMD) 70 
simulations are conducted under UHV and solvent-immersed electrode conditions, respec- 
tively. 

The majority of the calculations are performed using the Vienna Atomic Simulation Pack- 
age (VASP) version 4.6 7 - 1 "" 7 - 2 - and the PBE functional.— Modeling spinel Li x Mn 2 04 requires 
the DFT+U augmented treatment 7 " of Mn 3d orbitals. The U and J values depend on the 
orbital projection scheme and DFT+U implementation details; here U — J =4.85 eV is cho- 
sen in accordance with the literature.— A few PBE0 hybrid DFT functional calculations 
are also performed using VASP version 5. 2. 77 This more costly method is meant to provide 
spot checks of DFT+U results. A 400 eV planewave energy cutoff is imposed in all cases. 

Static geometry optimization and climbing image nudged elastic band (NEB) 78 barrier 
calculations are performed using a 10 -4 eV convergence criterion and 1x2x2 Brillouin 
zone sampling. The simulation cells are of 27.59 x 8.40 x 8.40 x A 3 dimensions and have a 
Li 6 Mn 2 o04o stochiometry, not counting the EC molecule. The bare lithiated spinel slab is 
generated as follows. A (100) Lii Mn 20 O40 surface slab is relaxed to its optimal geometry. 
The predicted 0.69 J/m 2 surface energy agree with the literature value.— To reach the target 
stochiometry, 4 Li atoms need to be removed. There are 210 possibilities, many of which 
are related by symmetry. The lowest energy configuration (Fig. [2]) is found to be the one 
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where a Li on each surface and half (two) of the Li atoms in the middle layer of the slab are 
removed. The lattice constants are kept at those of LiMn204. A larger system, LigMn2404s, 
is also discussed in appendix A. 

AIMD simulations feature 40x 11.88x 11.88A 3 simulation cells with a Li 12 Mn 4 oOgo ox- 
ide slab generated by expanded the oxide slab used in static calculations twofold. 32 EC 
molecules are confined between the two surfaces of the slab in the periodically replicated sim- 
ulation cell. T-point Brillouin zone sampling and a 10 -6 eV convergence criterion are applied 
at each Born-Oppenheimer time step. The trajectories are kept at an average temperature 
of T=450 K using Nose thermostats. Tritium masses on EC are substituted for protons to 
permit a time step of 1 fs. Under these conditions, the trajectories exhibit drifts of less than 
2 K/ps. EC configurations are pre-equilibrated using Monte Carlo simulations and simple 
molecular force fields, as described in an earlier work.™ We also report AIMD simulations 
of hydroxylated Mn0 2 surfaces in appendix B, using deuterium masses for H atoms and a 
0.5 fs time step there. 

AIMD potential-of-mean-force (PMF, denoted as W(R) hereafter) calculations apply 
composite reaction coordinates of the form R=Ri—R 2 . R covers the first two steps of the 
reaction. R± and R 2 are the distances between the Oe and Cc atoms and between Cc 
and Oi (on the oxide surface), respectively (Fig. [1]). In the next step, we apply another 
composite coordinate R' = R[ — R' 2 , where R[ and R' 2 are the distances between H and Ce, 
and between H and O2, respectively (Fig. [T]). Harmonic potentials of the form B Q (R — R Q ) 2 
and a series of 10 windows with a progression of R Q values span the reaction path R. B Q 
ranges from 2 to 4 eV/A 2 . Less windows are used along the R' coordinate. The statistical 
uncertainties in W(R) and W(R') are estimated by splitting the trajectory in each window 
into four, calculating the standard deviation, and propagating the noise across windows 
assuming gaussian statistics. 

To determine Mn charge states, we examine the electron spin on each Mn ion demarcated 
according to PAW orbital projections.— All calculations allow unconstrained spin polariza- 
tion. Bulk LiMn 2 04 crystal is antiferromagnetic (AFM). We have imposed an AFM ordering 
on alternate Mn planes along the (Oil) direction. The bare spinel surface slab is also AFM. 
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III. RESULTS 



A. Gas phase energetics 

The DFT/PBE functional predicts that gas phase disproportionation of ethylene carbon- 
ate into carbon dioxide and acetaldehyde, 

EC ->• C0 2 + CH3CHO, (1) 

is exothermic by 0.317 eV at T=0 K, not counting zero-point energy or finite temperature 
entropic effects which should further favor dissociation. This result is qualitatively consistent 
with a reported hybrid DFT predictions which however include entropy contributions. 79 We 
have not calculated the activation energy, but it must be high for metastable EC to be useful 
in battery operations. No electron is transferred in Eq. [TJ The corresponding dispropor- 
tionation products for PC would be CO2 and acetone ((CHs^CO). Therefore observation of 
CO2 and acetone/acetaldehyde in experiments may not always arise from electrolyte oxida- 
tion, as was already pointed out.— ^ They could be also be the product of EC/PC catalytic 
decomposition on battery material surfaces. 

B. Reaction intermediates 

Next we consider the intermediates associated with EC decomposition on Lio.6Mn204. 
Discussions of the activation barriers between them will be deferred to later sections. 

Figured (configuration "A") depicts an intact EC physisorbed on the (100) surface. The 
carbonyl oxygen (Oc) is the only atom strongly coordinated to the oxide, and is 2.16 A from 
one of the two surface Mn(IV) ions (among eight surface Mn, appendix A) in the simulation 
cell. The overall adsorption energy is 0.480 eV. This tilted EC geometry binds more strongly 
to the surface than one where EC is in upright. The overall magnetic polarization of the 
simulation cell is 2 fi-Q, If a net zero spin polarization constraint is imposed, the total energy 
increases by only 0.040 eV. This underscores the fact the energy differences associated with 
spin flipping tend to be small. 

In the first reaction intermediate B (Fig. [2b), the carbonyl carbon atom (Cc, Fig. CP is 
subjected to nucleophilic attack by an oxygen atom on the surface, becoming 4-coordinated 
and sp 3 hybridized. This bent EC species is reminscent of the bent EC - geometry during 
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method 


A 


B 


C 


D 


E 


DFT+U 
PBEO 


0.000 
0.000 


-0.011 
-0.107 


-0.100 
-0.260 


-2.057 
-1.911 


-2.144 
NA 



TABLE I: Energies of configurations A-D (Fig. [2]), in eV, computed at T=0 K using different 
methods. 

electrochemical reduction of EC at anode surfaces.— However, a maximally localized Wannier 
function 80 analysis performed on this configuration reveals that the EC carries no net charge. 
The configuration change is slightly exothermic (Tabled]). Note that two types of oxygen 
ions exist on the (100) surface even in the absence of adsorbed EC: those bonded to an Mn 
ion on the subsurface layer, and those that are not. The latter are more reactive, and the 
EC adsorption configuration is chosen so that the Cc atom can bind to the latter type. 

The next intermediate C (Fig. [2b ) features one broken Cq-Oe bond, with the Cc atom 
now reverting to sp 2 hybridization. Both the Oe atom on the broken bond and Cc become 
coordinated to the surface while the carbonyl oxygen (Oc) is detached. The surface Mn ion 
bound to the Oe atom now has 6-ionic or covalent bonds and transitions from a +3 to a 
+4 charge state, with an 3<i electron originally in its 3d shell migrating to a Mn(IV) on the 
surface (appendix A). 

In Fig. [2bl (product D), the proton on the Ce atom closest to the surface is transferred to 
a surface oxygen with a substantial gain in energy (Table [1]). Counting all Wannier orbitals 
centered closer to the EC C and O nuclei than Li .6Mn 2 O4 atoms as belonging to the EC 
fragment, EC now contains 32 valence e~ , implying that two of its e~ have been transferred 
to the oxide (appendix A). The oxidized EC fragment contains a +33|e| pseudovalent nuclear 
charge and a +|e| net charge. This oxidative step is apparently facilitated by the preceeding 
bond making/breaking events. 

The proton-accepting O ion on the oxide surface is coordinated to 3 other Mn ions. The 
average Mn-0 distance associated with this O atom becomes 2.13 A, an increase of 0.14 A 
from the value (1.99 A) prior to proton transfer. Thus EC oxidation on the spinel oxide 
surface has led to significant weakening of the Mn-0 ionic bond network on the Li2,Mn 2 04 
surface. We argue that this may contribute to Mn(II) dissolution^ (see also appendix B). 

Finally, as will be discussed in Sec. IIIIDl the configuration in Fig. Ebl spontaneously reor- 
ganizes during AIMD simulations of the liquid EC/Li .6Mn 2 O4 interface. An oxygen ion on 
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the oxide surface is pulled outwards, and the Oc atom re-coordinates to an adjacent surface 
Mn ion. Extracting the EC fragment configuration from an AIMD snapshot, transplanting 
it into the smaller UHV simulation cell to give a cell commensurate with those of Figs. |2^-d, 
and optimizing the geometry yield Fig. [2fe, with a further 0.087 eV gain in energy. Thus all 
steps of the reaction are exothermic or approximately thermoneutral (Table H]). The config- 
uration in Fig. [2^ maintains a +|e| charge in both the T=0 K optimized geometry and the 
AIMD snapshot, counting only the electrons on the original C and O atoms of the EC. 

DFT+U predictions depend on the value of U. This is especially a concern with the C— >D 
reaction, a redox process accompanied by Mn charge state changes. The redox potentials 
of transition metal ion centers have been shown to depend on DFT functionals in some 
cases.— - — Thus we have applied the widely used hybrid DFT functional PBEO, which carries 
no tunable parameters, to spot-check DFT+U results. The predicted exothermicities differ 
by less than 0.16 eV in all cases considered (Tabled]). In particular, the overall A— >D redox 
reaction at T=0 K differ only by 0.15 eV. The good agreement shows that the results 
are not strongly dependent on DFT functionals, and gives us confidence in DFT+U-based 
predictions for our system. 

C. Activation (free) energies: A to C 

Figure [3] depicts the zero temperature energy profile connecting A (Fig. [2h) and B 
(Fig. [2b), computed using the climbing image NEB method. The barrier is 0.240 eV. This 
small barrier is substantially less than the endothermicity associated with the EC bending 
motion on hydroxylated LiA102 surfaces,— presumably because EC binds more strongly to 
Li rr Mn 2 4 (100). 

The barrier between intermediates B and C has not been calculated at T=0 K. Instead, 
AIMD simulations of liquid EC/Li x M 2 n04 interfaces at finite temperature cover the A— >C 
steps. But first we discuss unconstrained AIMD simulations of EC liquid confined between 
Li x Mn 2 04 surfaces prior to any chemical reactions. Intact EC molecules have large dipole 
moments but are generally much more strongly coordinated to cations than anions. The 
8 surface Mn ions are observed to be coordinated to between 2 and 7 EC molecules, where 
we have used the arbitary but reasonable coordination criterion that Mn must be within 
2.5 A of one of the EC oxygen atoms. The mean value is 3.7 over the trajectory. As will 



9 



be shown, EC should decompose on a subsecond timescale on Li .6Mn 2 O4 (100). Therefore 
the distribution of intact liquid EC on the clean (100) surface is of academic interest.-^ The 
primary purpose of including the EC liquid is to provide a reasonable EC liquid solvation 
environment. 

PMF calculations are initiated from the unconstrained trajectory. We start with a 
snapshot where an intact surface EC molecule exhibits an adsorption geometry resembling 
Fig. [2^. As mentioned in Sec. [Til a composite reaction coordinate R=Ri-R 2 is applied, where 
Ri is the distance between Cc and Oe and R2 is that between Cc and Oi on the surface 
(Fig. [1]). If the simple coordinate R = R\ is used, the EC tends to desorb from the surface 
during the AIMD trajectory before bond breaking can occur. If R = R 2 , the reaction stops 
after configuration B is reached. 

Figure depicts the W(R) computed using this composite coordinate. The minimum 
associated with intermediate B is less stable than the reactant by 0.24±0.05 eV. This is less 
favorable than the —0.011 eV found at T=0 K (Table HJ. The barrier between A and B also 
increases to a still modest 0.48±0.05 eV. The reason is most likely the entropy change which 
is included in W(R) but not in Fig. [3^. The 4-coordinated intermedate B is evidently more 
constrained and entropically unfavorable than the loosely bound, flat EC molecule (Fig. [2^). 

The AIMD-predicted barrier between intermediates B and C is smaller than that be- 
tween A and B. The exothermicity between the W(R) minima associated with C and A is 
— 0.06±0.06 eV, which is similar to the energy difference computed at T=0 K (—0.10 eV). 
Thus, the main conclusion is that the liquid EC environment does not qualitatively modify 
the conclusions obtained under UHV conditions. 

D. Activation (free) energies: C to E 

Figure UK depicts the T=0 K energy profile between intermediates C and D. The reaction 
is exothermic by 1.957 eV and exhibits a 0.618 eV energy barrier. Ignoring entropic effects 
and assuming a typical molecular vibrational prefactor of 10 12 /s, such an enthalpic barrier 
permits a 18/s reaction rate. As a typical battery charging/ discharging rate is 1 C (i.e., takes 
about an hour), the predicted subsecond reaction rate means that EC oxidation proceeds 
readily during battery operations. 

Next we consider temperature/solvent effects. Figure Hb depicts the AIMD W(R') for 
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this segment of the reaction. We apply a second composite reaction coordinate R' = R^ — R^, 
where R[ and R' 2 are the distances between H and Ce, and between H and O2, respectively 
(Fig. [Th; Sec. [II])- The predicted liquid state W(R') barrier (0.54 eV) is slightly smaller 
than the gas phase value (Fig. Hk). The overall reaction exothermicity may be higher in 
the presence of liquid EC than under UHV conditions because the product is charged and 
stabilized by the high dielectric liquid environment. As the barrier is the main object of 
interest, we have not extended W(R') far into the product channel. It is sufficient to know 
that the reaction is strongly exothermic. 

Note that the electron transfer from C to D occurs right at the electrode-solvent interface. 
It is distinct from the long-range electron transfer mechanism experienced by "outershell" 
ionic complexes in classic electrochemical studies, and should consequently be much less 
influenced by the electric double layer. 

Configuration D spontaneously reorganizes to E (Fig. [2fe) after a 2 ps AIMD trajectory. 
We have not computed the free energy barrier between these configurations, but it must be 
on the order of k^T for the conformation change to occur on such a short timescale. 

E. Other products/intermediates are unfavorable 

Fig. [5k depicts an otherwise intact EC molecule with a proton transferred to the oxide 
surface. Compared to the intact adsorbed EC (Fig. [2^), this configuration is endothermic by 
0.398 eV. In contrast, removing a proton from the EC with a broken Cc-Oe bond (Fig. [2b) 
is exothermic by 2.057 eV overall. This shows that the non-oxidative, bond-breaking first 
steps strongly facilitate redox reactions on the spinel surface. 

Several possible decomposition products, devised from breaking a Ce-Oe bond at various 
stages of the reaction, lead to higher total energies and/or high barriers and are ruled out. In 
Fig. [5b, we attempt to break a Ce-Oe bond in the Fig. [2^ configuration to form a carbonate- 
like group in the first reaction step. The optimized metastable geometry is unfavorable by 
1.562 eV compared to the physisorbed, intact EC (Fig. [2^). 

Fig. [5b depicts a configuration resulting from cleaving a Ce-Oe bond in Fig. [2b without 
first transferring a proton to the surface. One product fragment, a acetaldehyde molecule 
which arises from internal proton transfer among the carbon atoms in the CH2CH2O frag- 
ment, is detached from the surface, while the other, a CO2 molecule, remains strongly 
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coordinated to the metal oxide. No electron is transferred. This configuration is exothermic 
by 0.646 eV compared to the intact physisorbed EC, but is much less stable than that in 
Fig. EH. (See also Sec. OTTftl ) We have not found a converged barrier smaller than 1.5 eV 
for this reaction. Hence the Fig. EH product channel is much more favorable. 

Figure |5H explores the possibility of Ce-Oe bond-breaking following proton transfer 
(Fig. EH)- The resulting product is far less stable than if the Ce-Oe bond remains intact 

(Fig. EH). 

In conclusion, we have not found a favorable reaction pathway that yields CO2 and/or 
acet aldehyde- like fragments. This suggests that the partially decomposed EC fragments 
chemisorbed on the (100) surface (Fig. EH &Efe) is stable under UHV conditions. Preliminary 
AIMD PMF simulations also show that breaking a EC bond in configuration D (Fig. EH) or 
configuration E (Fig. Efe) to Fig. [5H exhibits large barriers. 

However, the chemisorbed EC fragment (Fig. Efe) can potentially attack and react with 
other intact solvent molecules, potentially abstracting protons and/or electrons ("H~") from 
them. This has not been observed in our AIMD simulations which are of limited duration. 
In view of Refs. 



65 



and 



66 



it is also possible that further decomposition of the fragment can 



occur in the vicinity of anions in the electrolyte, which can act as proton shuttles. 



IV. DISCUSSIONS 



A. Relevance to battery operating conditions 

Our calculations focus on the clean Li r Mn204 (100) surface, which has been the focus of 
almost all theoretical studies.— 1 ^ Thin LI2CO3 films are known to exist on cathode oxide 
surfaces before cycling power.— >^ Since electrolyte decomposition products are observed on 
these oxides surfaces, it has been suggested that this film dissolves upon immersion in the 
electrolyte, during or prior to cycling power.— >2sL2£ 

The reaction products and pathways reported on (100) are not meant to quantitively 
reproduce all experimental spinel facets present. If the adsorbed, partially decomposed 
EC + fragments (Fig. EH, e) on Li .6Mn 2 O4 (100) are stable, they become part of the cathode 
surface film and may interfere with Li + intercalation. These products block the reactive 
sites, of surface concentration 1/17.6 A 2 , but do not form a continuous film to stop electron 
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tunneling from the cathode to the electrolyte. 

The fact that PFg is found to accelerate M112O4 degradation is not inconsistent with 
our findings. PFg may further react with, decompose, and removed the chemisorbed EC + 
fragments on the surface (Fig. [2fe) to form volatile organic products.— This may expose reac- 
tive sites on the oxide surface, allowing other EC molecules to react and potentially deposit 
more protons. This will be considered in the future. The H^O/PF^ pathway discussed in the 
literature may also occur simultaneously with our predicted EC H + -donation mechanism. 

As no counter electrode is present in the simulation cell, our AIMD simulations represent 
an open-circuit condition with a Lio.6Mn204 stochiometry which may be more related to bat- 
tery storage than charging experiments. For larger values of x, Mn(IV) has less tendency to 
reside on the oxide surface, and the reactivity of EC should decrease. This may be consistent 
with the observation that higher charging voltages yield more interfacial reactions.— 1 ^ 



B. Comparison with EC + /PC + decomposition predictions made in the absence of 
an electrode 



To some extent, our proposed EC decomposition mechanism (Fig. [2|i-e) is the reverse 
of that of Refs. and 1(36. Oxidation and proton removal from EC occur in the last, not 
first, step. A C-0 bond first ruptures on a EC initially physisorbed the Li a .Mn 2 04 surface, 
facilitating subsequent redox reactions. 



As in this work, Refs. 



65 and 



66l emphasize the role of proton removal from organic solvent 



molecules during oxidation. All these studies suggest that trace water may not be needed to 
initiate degradation of the cathode by creating acid molecules. The final products reported 



in Ref. 



65 



(which does not consider explicit electrodes) are coordinated to anions. In our 
work, anions are not present, and similar products are at best metastable. 

We predict preferential cleavage of the Cc-Oe bond, not Ce-Oe, in agreement with Xing et 
aVs results on the decomposition of isolated EC + and PC + .— ^ The final products predicted 
by Xing et al. are not observed because EC oxidation is accompanied by proton transfer to 
the surface herein. 
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C. More discussions on computational methods 

er transfer between the EC and Mn ions, and between Mn ions, are implicitly assumed to 
exhibit lower barriers than the rate-determining step in the bond-breaking processes. This 
follows from the DFT (and DFT+U) calculations, where e~ transfer events are implicitly 
assumed to adiabatically follow ionic trajectories. Non-adiabatic effects are not expected 
to dominate, because non-adiabatic polaron transport in spinel Mn 2 04 are not predicted to 
exhibit high barriers.™ The e~ transfer between configurations C (Fig. Eb) and D (Fig. [2]i) 
is unlikely to be rate limiting either, because the partially decomposed EC molecule is 
covalently bonded to the surface, yielding robust e~ conduction pathways. While hybrid 
functionals are generally more accurate than DFT+U for the C— >-D reaction which involves 
concerted electron and proton transfer, we have found that both methods predict similar 
energy changes (Tabled]). 

Zero point energy (ZPE) are not included in the energy changes reported in this work. 
Most bond-breaking events examined involve "heavy" second row atoms. In some cases, 
a proton is transferred from a carbon to an oxygen atom. As C-H and O-H vibrational 
frequencies are within 20% of each other, ZPE should not strongly modify the predicted 
exothermicities. Proton transfer barriers would be even lower if ZPE is included. 

Finally, we mention that the reaction pathway (Fig. [2]) is chosen through trial and error. 
A more systematic approach would be to apply metadynamics.— ~— The transition path 
sampling approach can also be combined with AIMD simulations to yield more rigorous free 
energy barriers." 

V. CONCLUSIONS 

This computational work focuses on the initial stages of ethylene carbonate (EC) de- 
composition on Li .6Mn 2 O2 (100) surfaces. In analogy with studies of water adsorption 
and dissociation on transition metal oxide surfaces, we perform calculations on clean, ideal 
surface models that can potentially be replicated in future experiments. Both ultra-high 
vacuum (UHV) and solvent-immersed electrode conditions are considered. Static (T=0 K) 
DFT calculations and finite temperature ah initio molecular dynamics (AIMD) simulations 
are applied under these conditions, respectively, yielding qualitatively similar conclusions. 
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The first step of EC decomposition on Li .6Mn 2 O4 involves Cc-Oe bond-cleavage but no 
oxidation. This slightly exothermic bond-breaking event enables a proton and two electrons 
to be transferred to the oxide surface in the next oxidative, strongly exothermic step. The 
activation energies for both steps are about 0.62 eV under UHV conditions and 0.54 eV at 
electrode-electrolyte interfaces, according to our DFT+U calculations. Thus oxidation of 
EC on the (100) surface is predicted to take place on a subsecond timescale. The Mn-0 
ionic bonding network is weakened by the EC proton transferred to the surface as well as the 
adsorbed EC fragment, suggesting that solvent decomposition and Mn and/or O dissolution 
from the oxide surface may be related. The fact that proton can come from organic solvents 
has been suggested by an experiment work using tetrahydrofuran.— Our work suggests a 
surface-assisted pathway for solvent oxidation and surface film formation at the open circuit 
Lio.6Mn204 voltage. No assumption or separate calculation of spatial inhomogeneity effect 
on the local dielectric constant is made in this work. 

The predicted EC fragment chemisorbed on the oxide (100) surface is likely not the final 
product. Our goal is not to reproduce the final distribution of organic decomposition prod- 
ucts observed in battery experiments. In particular, our calculations do not extend to the 
high voltage region which seems to be required for CO2, one of the proposed final products. 37 
Instead, this work represents the perhaps necessary first step of explicit modeling of elec- 
trolyte breakdown on oxide surfaces, and paves the way for theoretical and experimental 
investigations on other Li 3 ,Mn 2 04 crystal facets. 
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Appendix A: Charge states of Mn ions 

In Sec. IIII B\ changes in Mn charge states within the oxide slab are briefly discussed. 
They are explicitly depicted in this appendix. Fig. [H^ shows that 6 Mn(III) and 2 Mn(IV) 
ions in the Li 6 Mn 20 O40 slab reside on the (100) surface. The physisorption of an intact EC 
on a surface Mn(IV) ion does not change this charge state distribution (Fig. |6b). Upon 
breaking an EC C-0 bond, however, the O atom on that broken bond becomes strongly 
bound to a surface Mn ion, which loses an electron to the adjacent, formerly Mn(IV) ion 
on the surface (Fig. [Hb). The oxidative step transfers two e~ to the oxide slab. As a result, 
both a surface and an interior Mn(IV) turn into Mn(III) (Fig. EH). 

The possibility of an EC adsorbed on a surface Mn(III) ion is also considered. Upon 
optimizing the Kohn-Sham wavefunctions and the geometry, the Mn(III) coordinated to 
the EC carbonyl oxygen reverts to an Mn(IV), with the extra electron migrating to another 
surface Mn(IV), converting that Mn to a Mn(III) (not shown). This suggests that the charge 
states of surface Mn ions strongly depend on whether molecules are coordinated to them, 
at least on this (100) surface. 

We have also considered a larger, Li 8 Mn 24 048 slab, of the same lateral dimension as 
the Li 6 Mn 20 O40 model studied in most of this work. The tetrahedral site occupancies are 
determined by calculating the energies of all symmetry-distinct Li-configurations. This 
model features 8 Mn(III) ions, and in principle all 8 surface Mn could have been Mn(III). 
Instead, 2 of the 8 surface Mn remain Mn(IV), presumably because the surface layers are 
deficient in Li + . 

Appendix B: Hydroxyl and proton on (100) spinel Mn02 surface 

Interactions between H2O and Li x Mn204 surfaces are pertinent to how trace H2O in the 
atmosphere or electrolyte affects Mn 2 04 stability. (Water may also emerge from electrolyte 
decomposition.) Acids are known to degrade spinel oxide electrodes,- 1 ^ and water can be 
a source of H + . In this appendix, we consider the effect of decorating and functionalizing 
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Li x Mn204 (100) surfaces with H 2 0, H + , and OH~ groups. 

Figure [7^ depicts H 2 physisorption on the (100) surface at the x = 0.6 stochiometry. 
It exhibits a binding energy of 0.452 eV, comparable to EC physisorption. Breaking a O-H 
bond in water and displacing the proton to a neighboring surface O ion to form two new 
hydroxyl groups lead to a 1.027 eV binding energy after performing geometry optimization 
(Fig. [7b). Thus dissociative chemisorption is more favorable than physisorption by 0.575 eV. 
However, even this chemisorbed configuration may be metastable compared to one where 
the H + has penetrated into the subsurface layer (see below). 

Next, we consider the effect of multiple hydroxyl groups on the surface. We turn to 
Mn 2 04, without Li + intercalation, and decorate all surface O ions with H atoms and all 
surface Mn ions with OH groups. The geometry is relaxed to a local energy minimum 
(Fig. [7b). The first type of hydroxyl groups, originating from protonation of surface oxygen 
atoms, point outward or parallel to the oxide surface. The second type, hydroxyl groups 
formed by adding OH pendants to Mn, always point out of the plane. As there are more 
added H than OH, the charge-neutral slab now contains some Mn(III) ions. 

Starting from this optimized geometry, a 11 ps AIMD trajectory is conducted at T=450 K. 
The final configuration is shown in Fig. [7H. Some type 1 hydroxyl groups now point inwards, 
with their H + occupying sites that Li + normally reside in Li 2 .Mn 2 04. This is reasonable 
because the interior sites bracketed by 2 ~ ions should be electrostatically more favorable 
for cations and protons. In one instance, a type 1 hydroxyl group has abstracted a H + 
from a neigboring hydroxyl group to form a H 2 0. The oxygen atom of that water molecule 
coordinates to a surface Mn ion 2.3 A away in this snapshot, instead of the initial Mn-0 
ionic bonding distance (~ 1.9 A). That Mn ion now has only 4 ionic bonds in addition to 
coordination to the H 2 molecule, and should be more susceptable to chemical reactions 
and dissolution. These trends are consistent with the experimental observations that protons 
accelerate decomposition of spinel Li ;c Mn 2 04.- i ^ i ^i Despite the low temperature compared 
to experimental annealing conditions during synthesis,— evidence of H 2 loss is already 
observed. 

Type 2 hydroxyl groups appear stable over the course of the AIMD trajectory. As surface 
O atoms dissolves after H + attacks, however, even Mn ions attached to type 2 hydroxyls 
may become undercoordinated, exposed, and susceptable to attack by other H 2 moleules 
or the organic electrolyte. 
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FIG. 1: Ethylene carbonate (EC) molecule. The panel depicts an isolated EC weakly bound to 
the surface; three Li x Mn204 atoms are also depicted to explain the labeling scheme used in the 
text. Purple, dark blue, red, light blue, and white spheres represent Mn, Li, O, C, and H atoms, 
respectively. 
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FIG. 2: (a) Intact EC (A) on Lio. 6 Mn 2 04 (100) surface, (b) Intermediate B. Note that the EC 
Cc atom sits atop a surface oxygen ion that is not bonded to a Mn immediately below, (c) 
Intermedate C, with a broken Cc-Oe bond. The surface Mn(III) ion coordinated to the Oe now 
becomes a Mn(IV). (d) Intermediate D. A proton is transfer to a surface O ion. Two electrons 
are transferred; the surface Mn(IV) described in panel (c) reverts to a Mn(III). Configuration 
(d) reorganizes to (e) during AIMD simulations without further e~ transfer, (f) Most stable Li 
distribution in LieMn2o04o (100) surface slab, viewed along the (Oil) direction, about 3 unit cells 
abreast. 
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FIG. 3: (a) Energy profile linking configurations A and B, computed at T=0 K using the NEB 
method, (b) Potential-of-mean-force linking the A, B, and C free energy basins along the composite 
reaction coordinate R described in the text. 
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FIG. 4: (a) Energy profile linking configurations C and D, computed at T=0 K using the NEB 
method, (b) Potential-of-mean-force linking C and D along the composite reaction coordinate R' 
described in the text. In panel (b), the bottom of the well at R=l A has not been determined. 



26 




FIG. 5: Four decomposition products or intermediates not considered part of the main EC de- 
composition pathway, with zero temperature energies of +0.398 eV, +1.562 ev, —0.646 eV, and 
—0.376 eV relative to Fig. [2^, respectively. 
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FIG. 6: Magnetic moments of Mn ions corresponding to the bare Lio.6Mn204 slab, intact EC 
(Fig. EC with a broken C-0 bond (Fig. [2b), and oxidized EC (Fig. (2bl). Purple and dark-green 
spheres represent Mn(IV) and Mn(III), respectively. 
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FIG. 7: (a)&(b): Intact and dissociated H2O molecule on Mn02 surface. In panel (a), one of 
the water protons is obscured by the O atom. (c)&(d): The initial and final hydroxylated and 
protonated Mn02 (100) surface configurations along an AIMD trajectory, respectively. In panel 
(d), one OH group on the surface has abstracted a proton from another hydroxyl group and has 
turned into a H2O molecule; that oxygen is colored deep blue. 
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